Quantum mechanics in phase space: 
First order comparison between the Wigner and the Fermi function 
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The Fermi gp(x,p) function provides a phase space description of quantum mechanics conceptually 
different from that based on the the Wigner function W(x,p). In this paper, we show that for a 
peaked wave packet the gp{x,p) = curve approximately corresponds to a phase space contour level 
of the Wigner function and provides a satisfactory description of the wave packet's size and shape. 
Our results show that the Fermi function is an interesting tool to investigate quantum fluctuations 
in the semiclassical regime. 
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I. INTRODUCTION 

The Wigner phase space representation of quantum 
mechanics [1, 2, 3] is a very useful and enlightening ap- 
proach. It is of practical interest in the description of a 
broad range of physical phenomena, including quantum 
transport processes in quantum optics [4] and condensed 
matter [5], quantum chaos [6], quantum complexity [7], 
decoherence [8], quantum computation [9, 10], and quan- 
tum tomography [11]. Furthermore, the phase space ap- 
proach brings out most clearly the differences and sim- 
ilarities between classical and quantum mechanics and 
offers unique insights into the classical limit of quantum 
theory [12, 13, 14, 15, 16]. 

The Wigner phase space distribution function of a 
quantum state described by a state vector \yj) reads 

W(x,p) = ^J°° dye-fa^ (x + |) V (x - |) . 

(1) 

(For the sake of simplicity, we consider the case of a sin- 
gle particle moving along a straight line). The Wigner 
function provides a pictorial phase-space representation 
of the abstract notion of a quantum state and allows us 
to compute the quantum mechanical expectation values 
of observables in terms of phase space-averages. 

A different, almost unknown phase space approach is 
based on an old paper by Fermi [17]. As pointed out by 
Fermi, the state of a quantum system may be defined 
in two completely equivalent ways: by its wave func- 
tion ip(x) — (x\ip) or by measuring a physical quantity 
gF(x,p). Given the measurement outcome gp(x,p) = g, 
ip(x) is obtained as solution of the eigenvalue equation 
gF(x,p)vj(x) = gijj{x), where p = —ihd x . On the other 
hand, given the wave function ij){x) it is always possible 
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to find an operator g F {x,p) such that 

g F (x,p)^(x) = 0. (2) 
Using the polar decomposition 

ip(x) = R{x)e^ s{x \ (3) 

where R(x) and S(x) are real functions [R{x) > for 
any x], it is easy to check that identity (2) is fulfilled by 
taking 

9F (x, -ihd x ) = [~ihd x - s\x)] 2 + h 2 ^±. (4) 

R{x) 

Equation (2) implies that the corresponding physical 
quantity gp{x,p) takes the value g = 0. The equation 

gF ( x ,p) = \p-S'(x)f + h 2 ^- = (5) 

R{x) 

defines a curve in the two-dimensional phase space. In 
other words, as expected from Heisenberg uncertainty 
principle, we cannot identify a quantum particle by 
means of a phase-space point (x,p) but we need a curve, 
g F (x,p) — 0. Note that it is also possible to write equa- 
tion (5) in the form 



p ± =S'(x)±^j-n 2 ^-=mv M ±V^Q> ( 6 ) 

where m is the particle mass, vm = ^S' the Madelung's 

velocity [18], and Vq = the so-called quantum- 

mechanical potential [19]. Equation (6) locates two 
points, (x,p + ) and (x,p_), in the phase space for any 
x such that R"(x) < and R(x) ^ 0. 

The phase space Fermi function gp(x,p) and the 
Wigner function W(x,p) are at first sight unrelated. In 
particular, for the Fermi function there is no interpre- 
tation in terms of quasiprobabilities as for the Wigner 
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function. On the other hand, for a Gaussian wave packet 
the gp(x,p) — curve is an ellipse of area ttH [20, 21] and 
coincides with the phase-space contour level along which 
W{x,p) = Wmax/e, with W max equal to the maximum 
value of W. Different contour levels of W correspond to 
different "equipotential curves" gp = constant. The pur- 
pose of the present paper is to show that a similar rela- 
tion exists when the Gaussian shape of the wave packet 
is modified, provided the wave packet remains peaked. 
Finally, we will comment on the significance of our re- 
sults in the context of semiclassical approximations of 
quantum mechanics. 



II. GAUSSIAN PACKETS 



Let us first consider the Gaussian packet 



R(x) = G{x) = 



1 



S(x)=p x. (7) 



In this case Wigner function (1) reads 

W(x,p) = —e & , 
while the Fermi function is given by 

h 2 (x-x ) 2 2 h 2 

g F {x,p) = + {p-po) - ^2 ■ 

It is clear from Eqs. (8) and (9) that for Gaussian packet 
(7) we have 



(8) 



(9) 



1 

neh 



(10) 



Therefore, there is a one to one correspondence between 
the "equipotential curves" gp{x,p) = K and W(x,p) = 

1 — s2 K 

C, with C = ~^ . In particular, the gF — curve 

coincides with the curve W = = , with W max — 
W(x ,p ) maximum value of W. 



III. NON-GAUSSIAN PACKETS 

We now discuss the relation between the Wigner and 
the Fermi function when the wave packet is peaked but 
not Gaussian. Assuming a smooth, regular behavior of 
the packet around its maximum, we choose the analytic 
expression R(x) = CG(x)[l + P(x)], with P{x) polyno- 
mial [chosen so that R(x) > for any x] and C nor- 
malization constant which is irrelevant for our purposes, 
while S{x) is a polynomial. The Wigner function is then 
given by the Fourier transform (1) of 

V> (x + f ) r {x-*)= C 2 G (x+%)G(x- f) 

x [l + P ( x + %)+p( x -%)+p( x +%)p( x -%)] 



[S(x+S)-S(x-S)] 




FIG. 1: Plots of the Wigner function (left) and of the gF — 
curve (right, thick full curves) for R Gaussian and various S. 
Horizontal axis: —6 < x < 6, with x = i(x — xq). Vertical 
axis: —6 < p < 6, with p = f(p — po). From top to bottom: 

The thin dashed curves in 



\x 3 , and 



(11) 
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the right plots correspond to the contour level W = vy ™ ax of 
the Wigner function. Note that for the top plot such contour 
level exactly coincide with the g F = curve. 



A. Gaussian R 

We computed numerically W(x, p) for several functions 
P(x), S(x) and found that it strongly depends on S(x), 
while the dependence on P(x) is weak, as far as the wave 
packet remains peaked. Therefore, we first focus on the 
case R(x) = G(x). 
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Wigner functions for P(x) = 0, namely 
ip(x) = G(s)e* s < x >, 



(12) 



and several S(x) are shown in Fig. 1 (left plot) and com- 
pared with the corresponding gp = curves (right plots). 
Even though the wave packets in Fig. 1, with the ex- 
ception of the top plot [S(x) oc x 2 , corresponding to a 
squeezed state], are far from being Gaussian, the gF = 
curve still provides a rather satisfactory description of 
size and shape of the wave packet in phase space. 

This agreement can be explained by the following ar- 
gument. If we set P(x) = and consider the expansion 



S(x + *)-S(a 



= S'(x)y + ±S"i(x)y3 + 0(y 5 ) « S'(x)y, 
we obtain 

^(z+f)V* (x-^j^F(x,y)e^ 



(13) 



(14) 



where F(x,y) = G(x+%)G 



Therefore, the 



shift theorem of Fourier transform implies that, if 

F y [F(x,y)] = W G (x,p), (15) 

with T y Fourier transform with respect to the y- variable, 
then 



F(x,y)ei s '<- X M = W G [x,p - S'{x)}. 



(16) 



We can therefore conclude that the Wigner function cor- 
responding to wave vector (12) reads 



■5 2 [p-S'(x)] 2 
P 



W(x,p) « W G [x,p- S'(x)} = 

(17) 

Hence connection (10) between the Wigner and the Fermi 
functions approximately holds around the peak of the 
wave packet. The rather good agreement between the 
gF = curve and the contour level W = 1-v " ax is shown in 
the right plots of Fig. 1. We can conclude that, for peaked 
packets, the gF — curve is close to an equipotential 
curve of the Wigner function, enclosing a phase space 
area of the order of Planck's constant. 



B. Non-Gaussian R, S = 

We have seen numerically that the dependence of the 
Wigner function on P(x) is weak, as far as the wave 
packet remains peaked. As an example, we consider 
P(x) = l + a (x ~f o)2 , with a > so that R(x) > for any 
x, and 5 = 0. For a < \ the wave function tp(x) = R(x) 
has a single maximum at x = and the gF = curve 
again gives a good representation of the phase-space size 
and shape of the wave packet (see the top plots of Fig. 2). 
On the other hand, for a > \ the wave function exhibits 




CD 





FIG. 2: Plots of the Wigner function (left) and of the gF — 
curve (right, thick full curves) for P = 1 + ax 2 , a — 0.3 
(top) and a — 10 (bottom), and S = 0. Horizontal axis: 
— 4 < x < 4. Vertical axis: — 4 < p < 4. 



two maxima at x± = ±Sy2 — ^. For a 3> 1 the two 
maxima are well separated. Such a "cat state" exhibits 
non-classical features which impact on the structure of 
the Wigner function (see Fig. 2 bottom left, for a = 10). 
Since W(x,p — 0) is the autocorrelation function of R(x), 
then it reaches its maximum at x = 0. On the other 
hand, the marginal J dpW(x,p) = \tp(x)\ 2 = R 2 (x) ex- 
hibits a minimum at x — and this is possible thanks 
to the negative regions of W{x,p) (the white regions in 
Fig. 2 bottom left). For this cat state the gF = curve 
captures the two peaks at x± (see Fig. 2 bottom right) 
but not the non-classical phase-space structures of the 
Wigner function. 



IV. A DYNAMICAL EXAMPLE: THE QUARTIC 
OSCILLATOR 



The gF = curve is an interesting tool to investigate 
quantum fluctuations in the semiclassical region. As far 
as the wave packet remains peaked, that is, before the 
Ehrenfest time scale [22], its size and phase-space shape 
can be readily derived from the wave function (or from a 
semiclassical approximation of the wave function), with- 
out computing the whole Wigner function. Moreover, the 
distortion of the wave packet is directly related to S'(x), 
that is, to the Madelung's velocity. 

As a numerical illustration of the capability of the 
gF = curve to capture relevant features of quantum 
fluctuations, we follow in Fig. 3 the evolution of the 
Wigner and Fermi functions for the quartic oscillator 
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Refs. [23, 25] for the evolution of classical and quantum 
phase-space distributions. Details on the computation of 
the Fermi function are given in the Appendix. 

In Fig. 3 we compare the evolution the Wigner function 
with the evolution of the gF = Fermi function, for 
the quartic oscillator, starting from an initial Gaussian 
wave packet \a) (a\a) = a\a)). It is clear that, as far 
as the wave packet remains peaked, the gp — curve 
reproduces both size and shape of quantum fluctuations. 
This is the case for times smaller than the Ehrenfest time 
scale tE ~ [22], until which the centroid of the wave 
packet follows a classical trajectory. For longer times the 
Wigner function develops interference fringes, while the 
g F = function splits into several curves. We point out 
that the gp{x^p) = function (6) singles out only two p- 
valucs (or none) for any q. Therefore, it cannot reproduce 
the whole phase space structure of the wave packet when 
the Wigner function does not exhibit a single peak but 
a non-monotonous behavior along p [see the bottom plot 
of Fig. 3 (left)]. 




FIG. 3: Plots of the Wigner function (left) and of the qf — 
curve (right) for the quartic oscillator, with hX/uj — 0.01. 
Horizontal axis: — 10 < x < 10, with x = ^/^^-x. Vertical 

axis: — 10 < p < 10, with p = j^jP- From top to bottom: 
t/T = 0,0.4,0.8, 1.2, with T = 2-k/uj. The initial Gaussian 
distribution is centered in (x,p) = (3, 0). 

Hamiltonian 

H = hojcJa + h 2 \{a)Ya 2 . (18) 



Here, a = y h ^ (mcux + ip), u> is the frequency of the 
harmonic part of the oscillator and hX gives the strength 
of the nonlinearity. This model has been widely inves- 
tigated in the context of quantum to classical transi- 
tion [22, 23, 24, 25] and also used to explain important 
experimental results [26]. Model (18) is integrable, see 



V. CONCLUDING REMARKS 

In summary, we have shown that the phase space struc- 
ture of a peaked wave packet can be satisfactorily de- 
scribed by the g F = Fermi curve. In spite of the 
fact that the Wigner and the Fermi functions are at first 
sight two completely unrelated phase space descriptions 
of quantum mechanics, a link between them exists and 
is based on the shift theorem of Fourier transform. Such 
theorem also allows us to understand the shape of the 
Wigner function for perturbed Gaussian packets in terms 
of the Madelung's velocity. Our theoretical results, cor- 
roborated by numerical simulations for the quartic oscil- 
lator model, show that the Fermi function is an interest- 
ing tool to investigate the phase-space size and shape of 
quantum fluctuations in the semiclassical regime. 

While the Fermi function gp(x,p) fully determines the 
state of a quantum system, the extension of the results 
obtained in this paper to generic states encounters diffi- 
culties. Knowledge of the gp(x,p) = curve is in general 
not sufficient. The complete determination of the state of 
a system requires the extension of this curve to the com- 
plex p-plane. That is, consideration of the complex val- 
ues of p±, obtained from Eq. (6) when R" > 0, is needed. 

We then obtain S' = h 2 ^ = - ( £± ^)' \ from 

which the gp operator (4) and consequently the wave 
function ip(x) are determined. Any phase space descrip- 
tion of quantum mechanics necessarily involves features 
beyond classical intuition: negative quasiprobabilities in 
the case of the Wigner function, complex momenta for 
the Fermi curve. 
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APPENDIX A: FERMI FUNCTION FOR THE 
QUARTIC OSCILLATOR 

We consider as initial condition a Gaussian state (co- 
herent state for the harmonic part of the oscillator), 

I") = T,n=o c n\n), with c„ = e-sH 2 ^, a\a) = a\a), 
H\n) — E n \n), E n = hio + h 2 Xn(n - 1). The state of 
the quartic oscillator at time t is then given by \ip(t)) = 
S^Lo c ™ e ^ B "*l n )- ^ n * ne coordinate representation, 



1/4 



2 n/2 



<f> n (x) = {x\n) 



(Al) 



where H n denotes the n-th Hermitc polynomial. 

In order to compute the Fermi function, we write ip = 



ipu + iipi, so that 



R 



s' 



(A2) 



(A3) 



Finally, the gp — curve is obtained from (5). Note 
that in computing the derivatives of ipR and ipi, we took 
advantage of the relation H' n = 2nH n -i. This property 
of Hermite polynomials allowed us to avoid numerical 
errors in the computation of the derivatives R" and S'. 
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